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Abstract 

We provide an explicit analytical calculation that shows the asymptotic ap- 
proach of the one dimensional Caldeira-Leggett model to thermal equilibrium 
in the high temperature and weak coupling limit. We investigate a free par- 
ticle and a harmonic oscillator system, using both the Lindblad and the non- 
Lindblad type master equations for each case, and show that thermal equilib- 
rium is reached exactly for the free particle and aproximately for the harmonic 
oscillator, irrespective of the initial preparation of the system. We also gener- 
alize our calculation to higher dimensions. 

1 Introduction 

The Caldeira-Leggett model p] is a simple system-reservoir model that can explain 
the basic aspects of dissipation in solid state physics, and in the high temperature 
and weak coupling limit, can also account for quantum Brownian motion [2] [2] [I]- It 
consists of a particle, which is also called "the system", that interacts with a heat 
bath of simple harmonic oscillators through a linear term. 

In this study, we first derive the exact solution of the Caldeira-Leggett master 
equation in a novel way and then use our results to investigate the long time behaviour 
of the density matrix and the approach of the system to thermal equilibrium. 

We first briefly introduce the Caldeira-Leggett master equation in both the Lind- 
blad and the non-Lindblad form, and examine the stationary solutions of the Lindblad 
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form and their relationship to the thermal equilibrium density matrix (Sec. [2]). After 
introducing our basic mathematical tools, we explicitly solve the time evolution of 
the density matrix for a free particle system using the Lindblad form of the mas- 
ter equation, and investigate its behaviour in the long time limit (SecHJ). We then 
study the differences that would arise from using the non-Lindblad form of the master 
equation (Sec. [5]). Next, we repeat all the previous calculations for the harmonic oscil- 
lator system (Sec. [6]) and finally show that all our results are generalizable to higher 
dimensions (Sec. [7]). We end with a summary of our results and some comments 
(Sec. ED- 



2 Master Equation and the Stationary States 



Caldeira-Leggett model is defined by the following Hamiltonian 

Hr 



Hv H c H B 

_^ *> ^_ 



n " n v ' n 

where q, p are the position and the momentum operators of the system and {q n }, 
{p n } are the position and momentum operators of the bath oscillators, respectively. 
Hs is the system hamiltonian, H c is a counterterm for renormalization, Hi is the 
interaction term and Hb is the hamiltonian for the reservoir. For our study, the 
reservoir oscillators are described by an Ohmic spectral density with a Lorentz-Drude 
cutoff function 

J[OJ) = uo- 



7T n 2 + LU 2 ' W 

In the high temperature, weak coupling quantum Brownian motion case in which 
we are interested, this Hamiltonian leads to the following master equation for the 
reduced density matrix of the system [5] : 

j t Pt = -^[H s ,pt] - 2lW ^ 2 BT [q, [q,pt]] ~ i^[q, {p, Pt}\ ■ (3) 

In order for this time evolution to be a quantum dynamical semigroup, one can add 
a "minimally invasive" term, — 8m 1 BT \p, \p,Pt]}, which is negligible compared to the 
other terms in the high temperature limit and which brings the equation into the 
Lindblad form [5] 

J t Pt = —[HsiPt] ~ 2j ^h hPt]] - t^ph M) ~ {P,Pt}} (4) 
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with the relaxation rate 7 and the Lindblad operator 



which gives the Lindblad evolution 

^Pt = ~[H s + ^(qp + pq),p t }+j\Ap t A^ -^Ap t -^p t A^A^j . (6) 

Certain conditions have to be met for ([3]) and (jl]) to be valid [5]: 

1. The typical time scale over which the state of the system changes appreciably, 
Tr ~ I/7, should be much larger than the typical decay time for the correlation 
functions of the bath oscillators, t b , which translates into: 

hrj < min{/ifi, 2-nk B T}, (7) 

The origin of this condition is related to the approximations that make the 
master equation Markovian and its details can be found in [5J. 

2. The typical system evolution time, 1/ujs, should be large compared to t b 

hw s < min{frft, 27ik B T}, (8) 

This condition is required for the validity of the introduction of the "minimally 
invasive term" as well as for other steps in the derivation of the master equation. 
As we expect p ~ muosq for the typical momentum and position values, the ratio 
of the momentum double commutator to the position double commutator in (J4]) 
is at the order (fiius/k B T) 2 , thus it vanishes under this condition, making the 
Lindblad form of the master equation valid. 

The first step in understanding the implications of the master equation is studying 
the stationary solution, i.e. the solution with dpt/dt = 0. The general expectation 
is that, in the long time limit, the density matrix is going to reach this solution 
irrespective of the initial condition. For the non-Lindblad case, ([3]), and for a poten- 
tial V(q) whose spatial variations are small, the stationary solution in the position 
representation is 

/ I I \ ~ at ( ^((gi + 92)/2) mk B T{ qi -q 2 f \ 
< gi |p|g 2 >-iVexp( — — ) , (9) 
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where N is a normalization constant [5]. Moreover, for the case of the quadratic 
potential or the free particle, this equation is exact. 

For the free particle, (V(q) = 0), gives the exact thermal equilibrium state, 
which can be obtained transforming the familiar expression 



p 2 

(pi\Pth\P2) = (pi\Ne 2mk BT\ p , 




Pi = P2 
otherwise 



(10) 



into the position representation, using the Fourier transform. 

For the harmonic oscillator, V(q) 
matrix is given by [6] 



\muj 2 q 2 , the thermal equilibrium density 



(Qi\p\Q2) = N exp 



2h tanh (hu / k B T) 



(<l 2 i + <&) 



muj 



h sinh (ftw/k B T) 



Qi Qs 



(11) 



which clearly does not agree with (Q. The source of this diagreement can be traced 
by expanding the exponent in ffTTl) in powers of frw/k B T. To the leading order 



(<7i|p|<?2) = N exp 
= N exp 



mk B T ( \ 2 (qi + q-i 



2h 2 



~8k B T 



k B T 



(9i + ^ 



4 

mk B T 



mk B T 
2h 2 



1 4~ 



2h 2 



Q2j 



exp 




which agrees with (Q on the diagonal (qi = 52), ignoring the higher order terms, but 
does not agree with it if we move away from the diagonal where q\ — q2 is comparable 
to qi + (/2 in magnitude. This should not be surprising, since the term in the exponent 



that causes the difference, 



24fc s T 



(<?i — <?2) 2 , is at the order {ftw /k B T) 2 compared to 
2 ^"(?i _( ?2) 2 - Also, we expect mu 2 q 2 ~ fku, which means (gi -g 2 ) 2 ~ -j^f- We 
already stated that our master equations is valid in the regime where hw/k B T 1, 
thus it is expected that such terms are not captured by the stationary state of the non- 
Lindblad master equation. In other words, by agreeing to use (EJ) or (0J to investigate 
the Caldeira-Leggett model, we give up the prospect to have sensitivity to any higher 
order terms in hu/k B T 

To summarize, the stationary solution of the free particle Caldeira-Leggett master 
equation in the non-Lindblad case is the exact thermal equilibrium density matrix 
for the free particle. Thus, reaching the stationary state and reaching the thermal 
equilibrium are equivalent for this case. However, for the harmonic oscillator, the 
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stationary solution is not the thermal equilibrium, but rather, the two agree on the 
diagonal in the high temperature limit. High temperature limit is essential to the 
derivation of the master equation itself, thus, we can say that the stationary solution 
captures the thermal equilibrium as best as possible for the Caldeira-Leggett master 
equation. 

Our previous discussions suggest that the "minimally invasive term" should not 
affect the physics significantly. We are not going to explicitly calculate the stationary 
solutions with this term, but our subsequent analysis will make it clear that the long 
time behaviour is identical to the non-Lindblad case. 



Different techniques have been used to solve the Caldeira-Leggett model, especially for 
the harmonic oscillator case [7] [8] [9]. Among these, [9] is the closest to our approach, 
but it only discusses the non-Lindblad form of the master equation and cannot be 
readily generalized to the Lindblad master equation (see appendix[D]for more discus- 
sions of this). 

For our analysis, we employed the techniques given in the appendices of [10] 
and [11], where they are used for similar but different purposes. The main feature of 
our technique will be the use of the Wigner function and the Gaussian ansatz. Details 
of the following calculations can be found in appendix lAl 

The characteristic function associated to the Wigner function is defined as 



where {\y)} are the momentum and position bases, respectively. The inversion 

formulae in the momentum and position basis are 



3 The Green's Function 




(13) 




(14) 
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respectively. After some algebraic manipulations, the time evolution for p t (k,x) for 
the case of V(q) = is given by 

d _ . . 1 , d _ . . 2^mk B T 2 7 2 _ 9 „ , . 

di PtiK X) = m k Tx Pt ^ X) ~ X X) ~ 8^f k Pt{K X) ~ 27 X Yx Pt ^ X) 

(15) 

We will employ the Green's function method to solve (115p 

p t (k,x) = / dk dx G(k,x,t; k ,x 0l O)po(k 0j x ), (16) 



where G(k, x, t; ko, Xo, 0) is the Green's function which is defined as the solution of 
( TT5l) that satisfies the initial condition 

\im.G(k, x, t; k , xq, 0) = 6(k — k )5(x — x ). (17) 



The key observation is that, under ([15]) . initially gaussian states remain gaussian. 
If we make the ansatz 

Pt(k, x) = exp{— c\k 2 — c 2 kx — c 3 x 2 — ic^k — ic 5 x — c 6 }, (18) 

then the master equation leads to 

— 6\k 2 — c 2 kx — c 3 x 2 — ic A k — ic 5 x — c' 6 

1 2'jmkBT 2 7 2 

= —k(-c 2 k - 2c 3 x - ic 5 ) — x - - — — -k - 2^x(-c 2 k - 2c 3 x - ic 5 ) . 

m nr ^rnksT 

(19) 

Upon equating the coefficients of the independent terms, we have the following system 
of ordinary differential equations: 

°2(t) , 7 



m SmkBi 

C2{t) = 27c 2 (t) 

m 

. 2^mk B T 

c 3 (t) = - 2 4 7 c 3 (t) (20) 



c 5 {t) 



m 

c 5 (t) = -2 7 c 5 (t) 
c 6 (t) = , 
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which are readily solved using basic techniques: 



ci(t) = Cl (0)+c 2 (0)— - + c 3 (0) 



2m7 

C2 (t) = C2 (0)e- 2 ^ + c 3 (0) 



4m 2 7 2 Sh 2 m^ 2 



(r t 2 + 2r 



k B T 7 



mksT 



2h 2 

c 4 (t) = c 4 (0)+c 5 (0) 1 ' 



(c 3 (0) 

r, 



r t e- 2 ^ | fc B r. 

m7 2/i 2 7 
mksT 
~ 2h 2 



-47* 



2m7 



c 5 (t) = c 5 (0)e- 2 ^ 
c 6 (t) = c 6 (0), 

where I\ = 1 — e~ 2lt . 

An initial gaussian of the form 

p k o OXO ' en (k,x) 



1 



-i(fe-fco) 2 e -^(^o) 2 



(21) 



(22) 



has the limit 



x) > 6{k _ ko)s{x _ Xq) 



and evolves in time to 



(23) 



/ fc B T i 7 \ l,2-> S- 

x e «b t / e 



(24) 



We find the Green's function given in j3] by taking the e — > 0, rj — > limit: 



G(/c, x, £; fco, xo, 0) = 5(A;q — fc)5 ( xo — (xe~ 27 * H ) 

2rrry 



X g \2h 2 m~/ 8mk B T 



\k 2 t 



fc2 r t +2r t _ fc £j_ 2 (1 _ e -4 7 t% 



(25) 
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4 Time Evolution of the Free Particle System 

The previous section set up all the background we need to solve for the time depen- 
dence of the density matrix. We first insert ( |25i) into ( {TBI to find pt(k,x) 



p t (k,x) 



g \ 2fi 2 m7 SmkgT 



fc2 l 1 i_ f lt_ fca ,_ ! _ L _ 3 , 2(1 _ e _ 47t) 

4m 2 j z m~l <• ' 



p {k,xe- 2 ^+ V,k 



2mj 



(26) 

In the long time limit, we expect to find that the density matrix will evolve 
to the stationary solution, which, as we have shown, is diagonal in the momentum 
representation. So, we will solve the time evolution of the density matrix in this 
representation: 



2nh 
1 



dx e ^ 2 >p t (jp 2 — 



Pi) x ) 



g V 2ft 2 7717 SmkgT 



(P2~pift 



2 (r^+2rj)fc B T 



dx e * ( 2 j e 



. r t (P2-Pl) 3 ._ (1 _ e -4-yt )a .2 



i?(t) + 



27ih 

x J dx e 
= R(t) + 



Po(P2 - Pi, xe 7 + 

r t (p 2 -pi)x 



2 7 t , r t (p2-p!) 



2m7 



Po(P2 "Pi, 



2m7 



-47*-, 



) mk B T ^.2 / i(Pl+P2) | r t fc B T 



2ft 1 2/v 



+-£d=-(»-Pi) * 



2wmk B T y/i - e" 4 ^ 



V2ft^m 7 + 8'«fci 3 Ty l - p 2 Pi J I gft 2 m7 2 U> 2 Plj 



x e 



2(l-e- 4 7*)mfc s T 



/ '(P1+P2) I T t k B T \Y 
\ 2fi + ^X^^ 2 PlJJ 



Po(P2 "Pi 



r f (p2 - gi) 

2m7 
(27) 



where 



i?(t) 



g V 2ft 2 m7 SmkgT 



(P2-P1) 2 * J r % + »f^ BT ^-Pl) 2 



2vrft 
x l dx { e 



(3 Sh-^rwy 



(i-e-^W ^ _ / M + r|^ (p2 _ pi) )a 



~ / -2 7 t , r *(P2 "Pi; 

Po(P2 - Pi, are 7 + 



2m7 



\ ~ / r t (p2 - pi)/ 

) - Po(P2 - Pi, ^ ) 



2m7 



(2? 
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As long as po(k, x) is bounded and well behaved around (p 2 — Pi, |~ ), R(t) vanishes 
in the long time limit. 

Then, the main observation is that as t — > oo, (pi\pt\p2) vanishes except for the 



case pi — P\ = 0, due to the term exp 
fact that 



k B T 



+ 



2h 2 m-y ~ 8mk B T 

p t (0,0) = tr p t = 1 , 



(p 2 - Pi) 2 ^ 



we finally reach 



(Pl|Poo|P2) 



2irmk B T 





e 2mfc s T p x = p 2 

otherwise 



. Using the 
(29) 

(30) 



which is the stationary solution of the master equation and the thermal equilibrium 
density matrix of a free particle (see Sec. (j2J)). Thus, the thermal equilibrium is indeed 
reached in the long time limit for the free particle in the Caldeira-Leggett model. 

All of the asymptotic behaviour of the density matrix can be read off from (I2T|) . 
Let us start by analyzing the remainder term R(t). For large times and for the values 
of x where the integrand is significantly different from 0, we can series expand po 
through the leading term 



PoG°2 - Pi, are 2lt + 



Ft(p2-Pi, 
2m7 



Po(P2 - Pi 



r i (p 2 -pi)x dp (k,x) 



2 7717 



dx 



(31) 



xe 



-27* 



given that po is well-behaved around x = (p 2 —p\)/2m^ and dpo(k,x)/dx does not 
vanish at the point (p2 — Pi, ) ■ Under these conditions, we can calculate R(t): 

-37* 9pp{k,x) 



R(t -> oo) 



g V2ft 2 m7 SmkgT y 

2irh 



<9x 



(p 2 - P1 ,to)) 



x 



dec 30 c 



(l-.-^}mi B T j 
2^ S " 



2h +^r£-(P2-Pi) 



^(Pi+pO , T\k B T \ dp (k,x) 

+ ~ tttt; — (P2 -Pi 



^2vr((l - e-^)mk B Tf \ 2H 2ft 2 7 ^ ^7 9x 



x e 



2(l-e- 4 T')mfe s T 1 



x e 



-(5^ + 8^fer)(P2-Pi) 2 t 27t 



/(Pi,P2) e v5fe + ^feW (p2 ~ Pl) * e' 2 ^ 



(32) 
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where 

,/ \ ti 1 fi(pi+P2) , k B T A dp (k,x) 

f{PuP2) = , , = = + 7TT^(P2 - Pi) 



yJlTs{mk B TY V 2/1 2^ 7 V " " 7 ^ 



R 

x e 2(l-e-W)mk B T\ 2ft ' 2ft^ 7 V- ^ " gfi 2 m7 2" ^33^) 

This means, for the non-diagonal elements, R(t) is dying much faster than the non- 
R(t) term in (1271) . due to the extra exponential factor of e -2 " 7 *, and is negligible. Then, 
the long time behaviour becomes 



1 {V2-V\) 2 k B T (Pi+P 2 ) 2 '(P2-P1) _ P2 ~ Pi 

(pi\pt^oo\P2) = \ — — e 2ft2 ™^ e 8m 's T e P0G02-P1,— ) 

V 2nmk B T zmj 

X e -(2^+8^T)fe-Pl) 2 i ^ ^ ^ ^ 

which is manifestly exponential. The time constant for relaxation is 

_ ( 7 A; g T V 1 1 _ 1 2ft 2 m 7 

P1 ' P2 \8mk B T ^ 2ft2 m7> / (p 2 - Pl ) 2 1 , f fry V k B T{p 2 - Pl f 1 J 

which leads to the following observations: 

1. In the momentum representation, relaxation to the thermal state becomes faster 
as one moves away from the diagonal in the density matrix. 

2. At a first look, relaxation is faster in both high and low temperature limits 
and is slowest at T = H r y/2k B , but remembering that Hj/k B T <C 1 should 
hold for the validity of our master equation, only the high temperature limit is 
meaningful and the ^j%mk B T term can be neglected: 

2m/i 2 7 2ft 2 

k B T(p 2 - Pi y D( P2 - Pl y [M) 



Here, D = k B T/m~y is the diffusion constant, which can be readily demonstrated 
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p=q 



using (J2 

(q 2 ) = tv(q 2 p t ) 

dp (p\q 2 p t \p) 

= -h 2 J dp (J^(q\p t \p) 

= — h 2 [ dpi —2 ( — ^ h - — 7 ml t (q\Pt\p) + time independent terms] 

J \ \ 2/7/7777 8mk B l J J 

« (t-Kx), (37) 

7717 

where we again used the fact that the trace of the reduced density matrix is unity. 
This is the behaviour of the classical Brownian particle, and shows the connection 
between the decay constants and the diffusion coefficient. 

For the diagonal elements, i.e. p\ = P2, there are two sources of correction to (13Qj) . 
one arising from R(t), and the other from the corrections to the non-i?(t) term due 
to the finiteness of t. The latter is at the order e -47 *, as can be seen from (127)) . so the 
leading correction comes from R(t) 



e~ 2 ^ ] , (3* 

(0,0) 



2nmkBT I mksT dx 

which shows that the relaxation time for the system is indeed tr ~ 1/ 7. 

5 The Time Evolution for the Non-Lindblad Mas- 
ter Equation 

The originally derived master equation for the Caldeira-Leggett model, ([3]), was not 
in the Lindblad form, lacking the "minimally invasive term" — 8tt J bT [P; [p>P*]]- The 
only difference not having this term makes in our previous analysis is that, when 
solving the differential equations for the coefficients of our Gaussian ansatz, the first 
equation reads 

m = C M , 09) 

777 

which is missing the ^/RmkBT term. All other differential equations for c;(i) are the 

7(P2 —Pi) * 

same as before. This in turn leads to the formally simple change that e 8mfc s term 
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is not present in (J2 
(Pl\Pt\P2) = R(t) 



I 1 fe B T c (r?+2r t )fe B T, , 2 

w 2nmk B T y/i - e -*fi y ' 



x e ^^v a£ ^ 2+ ^'— ■>)■ _ £ £ fe-p 1 ) 

2m7 

(41) 

This means: 

1. Any initial density matrix still evolves into the form in (1301) in the long time 
limit. So in the long time limit, the non-Lindblad and Lindblad equations give 
the same result, which is thermal equilibrium. 

2. Non-diagonal terms still vanish in the long time limit due to the e m2 ™t (P2 Pl> 
term, but the decay at first seems to be slower and the time constant is 

r - 2m ^ 7 (42) 

as opposed to (135|) . Nevertheless, remembering the condition thy <C min{M7, luk^ 
we can see that the missing term in the non-Lindblad case is negligible, so the 
change in the time evolution for the non-diagonal terms is negligible. 

7(P2~Pl) 2 i 

3. Since the e 8mfc s term was constant and equal to 1 for the diagonal matrix 
elements in the Lindblad form, not having this term does not have any effect. 
The time evolution of the diagonal elements of the density matrix is exactly the 
same as before. 



6 Time evolution of the Simple Harmonic Oscilla- 
tor System 

For a general potential 



V(q) = Y i a m q m , (43) 
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we cannot follow our method that solved the free particle case, since the equation 
contains terms of order higher than 2, and in that case Gaussian solutions are not 
preserved. Still, it is possible to have an analytical solution for the exceptional, but 
important, case of the harmonic oscillator. For a system under the potential 



V(q) = -mtu 2 q 2 



(44) 



all the calculational steps are very similar to, but algebraically more complicated 
than, those for the free particle, and they lead to (see appendix [B]) 



G(k, x, t; k , x , 0) 
e -(7-n)* A t 



6 k 



2/i 

x exp 



((yU coth fit + 7)& — muj 2 x) ] 5 ( xq 
k B T 



2h 2 muj 2 



Mi(t)k 2 + — — M 2 (t)kx + m 2 uj 2 M 3 (t)x 

7 



(fi coth fit — rfx + 
> (45) 



and 



Po 



Pt(k,x) = exp 

' e -(7-M)t A t 

2f~i 



mk B TM 3 (t) 
2h 2 



((yU coth fit + j)k — muj 2 x) 





k B T 


| expj 


2h 2 muj 2 



M x {t)k 2 + M 2 (t)kx 



2/x 



7 



(/x coth /it — 7)2; + 



fc 



(46) 



where /i = \f^ 2 — uj 2 , A t = 1 — e 2>lt and Mj are dimensionless functions with the 
asymptotic behaviours 



Mx(t -> oo) = 1 + 



fry 
2k B T 



+ 



M 2 {t -> oo) 



Af 3 (t -> oo) = 1 



1 / fay 



8 \ k B T 



4k^T 



4k B T 



+ 0[e 



-Re{2~/-n}t 



+ 0(e 



-i^e{27-At}^^ 



(47) 



For the case of the harmonic oscillator, we will use the position representation where 
the thermal equilibrium density matrix is given by (fTTl) . Using the inversion for- 
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inula (fT4|) . we have 

(QiIpM = 7TT / dke-^ k ( q - 1 ^)~p t (k,q l - q 2 ) 
27Tn J 

1 mk B TM 3 (t), s 2 f kgTM^t) u2 ( i(q 1 +q 2 ) , fcfl ™2 (fHgl Zgg) ' 

= ifo (t) + — e ^ — (9l - 92) / dke~ 2 ^ 2 - 2 -{^^ + ^ , 



with 



(48) 



RHo{t) = Y~fr e dke 2fe2 ™- 2 v 2ft / (p (k,x) 



where 



e -(7-/*)* A 4 

fc' = : ((/x coth /it + 7) A; — mu; (gi — g 2 )) 



2/j 

e -(7-M)* A, 



2/i 



' (/icoth/it -7)(gi - g 2 ) + — 1 • 



rn 



and we again used the fact that p t (0,0) = 1. 

Assuming that po is well behaved around (0,0), 



lim R HO (t) = . 



t— >oo 



(50) 



This means, in the long time limit 

(ffilPoolga) 



^ 2ttA; b T 



fry 



2ksT 



ftjjj 
4^T 



i -1 



x e 



,4k„T \4k n T \2k„T , 



2k p T ^\4k„T, 



1+ (^) 5 



x e 



16fc H T U H T, 



/17 



, 2fc R T I 1 V 4fc R T / 



(51) 



This result does not exactly agree with the thermal equilibrium density matrix (fill) 
or with the stationary solution to the non-Lindblad master equation (Q. Neverthe- 
less, conditions for the validity of the master equation imply that hj/ksT <C 1 and 
hw/ksT <C 1. Also, mu 2 q 2 ~ hw and 
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for the typical length scales we encounter in the harmonic oscillator system. Putting 
all these conditions together: 



<*W*> - t/^ e -iS£c« + *>V^<».->' . (53) 

So, the density matrix approaches the stationary solution (j9j) approximately, which is 
the expected result since was derived for the non-Lindblad master equation. Also, 
remember that the stationary solution agrees with the thermal equilibrium density 
matrix of the harmonic oscillator only on the diagonal and to the leading order in 
hu/k B T. 

Let us now analyze the corrections to the infinite time matrix elements. As- 
suming p is well behaved around (0,0) and defining dp (k,x)/dk\^ 00 ^ = D k and 
dpo(k, x)/dx\, Q ) = D x , the remainder term has the long time behaviour 

1 _ ^ B TM 3 (t) {qi _ q2)2 e ~(7-M)t At f f ^ _ ^ Dx 




R HO {t)^—e sp— 1 [(jjLcothtit + j)D k , 

_ k B T Ml ( t) k2 _ t ifa+gg) k B TM 2 (t)( n -gg) y / -muj 2 D k + (ucothut - 7)^ . 

Gt/Ce 2h2 miJ 2 V ft 2 7 / fc -| (fli — Oo 

^ (^coth/ii + ^Ac + ^f 

my 2 mk B TM 3 (t) , , 2 fe 2 m^ 2 / ifa+gj) | k B™2 CO (91 ~ 92 ) \ 2 
g 57T^ Wl-92; g2fc s TM 1 (i)V 2h 1 7 

27r/c s TMi(t) 
x ((pcoth/it + 7 )Ac + (D x /m)) 

/ hmw 2 (Aik B T(qi + q 2 ) - h^/(qi - q 2 )) -mu 2 D k + (/icoth/it - j)D x . 
X V 8^T2 + (^ol^+7)^T(^7H ~ 

x e- (7 -^* (54) 

In the long time limit, again there are two sources of correction to the infinite time 
values of the matrix elements. The first is from Rho which is exponentially small 
compared to the matrix element at the infinite time limit by e - & {7-M}i ; as seen 
in (I54p . The second correction comes from the corrections to the non-Rno term due 
to the finiteness of t in (l4~8l) . which mainly arises from Mj(t) — Mj(oo), and is at 
the order e - ^ 27- ^* as seen j n ( |47|) . Thus, after a sufficiently long time, the latter 
correction becomes negligible compared to the former unless Rho vanishes, which is 
only possible for a special combination of values of D k and D x . That is, the leading 
correction to flBTl) is (1541) . which dies out with a time constant of l/i?e{7 — /i} 

If we were to use the non-Lindblad master equation, than the only change would 
be that we would have no ^ /RinksT term for i\ in (!73j) . Effects of this can be easily 
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traced by having fcgT — > oo and H — > oo while keeping f? /ksT constant. This way, 
7/8mfcsT — > 0, with all other coefficients in the differential equations remaining the 
same, giving us the non-Lindblad equation. Then 



fojj h oo 



ksT ksT h 

hy h 2 7 







ksi Kb J- h 

In this limit, terms of 0{h^ /ksT) or O(hwfkBT) in the expressions for Mi, ( 1471) . 
vanish. Thus, without any approximations 



MPM = e -^<«.+»>V^<«--^ (56) 

as expected. This result was also reached in [9] , where they also show that the density 
matrix is diagonal in the energy basis of the system harmonic oscillator, supporting 
the idea of the pointer states of Paz and Zurek [12] . 

For the non-Lindblad master equation, The leading correction to the thermal equi- 
librium density matrix at large times is the same as the Lindblad case and originates 
from R H o(t). Thus it is 0(e~ Re{ ^~^ }t ). In short, non-Lindblad time evolution is 
equivalent to the Lindblad case within the sensitivity of our calculations. 

Note the following observations for the harmonic oscillator system 

1. If uj > 7, fj, is imaginary, and the relaxation time constant is l/i?e{7 — //} = 1/7. 

2. If uj < 7, ji is real and fi < 7, so relaxation still occurs but at the slower rate of 
e -ii~\/i 2 -u 2 )t _ The ra t e decreases as - approaches 0. 



7 Caldeira-Leggett Model in Higher Dimensions 

For a quadratic potential and a generalized dot product interaction term of — Ylt=i 1i K n,iQn,i 
the Caldeira-Leggett Hamiltonian for a free particle or a harmonic oscillator in d- 
dimensions is separable 

H = lt{^ + + «? E ^jr i + E (ik<i + ~ 

A II on fb*b an \ / 
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This means, we can write p t as a linear superposition of terms like <8> ...pf 1 ... <g> 
where 

^p? = -^[HsfrP®] - 2lim ^ BT [qu h,pt ] }} - 8m fe T \Pi, [Php?]} - ijbh {pupP}} ■ 

(58) 

We can solve for each dimension separately using its own parameter ji, take the tensor 
product and finally superpose such solutions to find the overall density matrix. So 
our results for the one dimensional case readily generalize to any higher dimensions. 



8 Conclusions 

We have solved the complete time evolution of the Caldeira-Leggett model for the 
free particle and the harmonic oscillator, and showed that the reduced density matrix 
of the system approaches the exact thermal equilibrium for the free particle and an 
approximate thermal equilibrium for the harmonic oscillator in the long time limit. 
Moreover, we studied the leading corrections to the infinite time behaviour of the 
matrix elements of the density matrix at finite times, and calculated the time scales 
for the decay of these corrections. We discussed the deviations from the thermal 
equilibrium for the harmonic oscillator in the infinite time limit and showed that 
they follow from the approximations we introduced during the derivation of the master 
equations. 

Our calculations explicitly showed that the typical time scale for reaching the 
thermal equilibrium is ~ I/7. Also we demonstrated that, within the limits of our 
master equations (J3j) and (Tj0), the Lindblad and non-Lindblad forms give the same 
time evolution. 

We note that the methods we used, the characteristic function of the Wigner 
function and the Gaussian ansatz, can be generalized to any equation of type ( JTBl as 
long as they do not lead to nonlinear variations of ( 120|) (see appendix [A]). We gave a 
treatment of the most general equation in the appendices. 
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A The Master Equation for the Characteristic Func- 
tion associated with the Wigner Function for 
the Free Particle Case 

In the position representation, #2) = (lilPtlte) , Q is 

dp t (qi,q2) 



Of 



ih ( d 2 d 2 \ i . t r / ^ 2 im k B T 



h 2 -i ( 8 d \ 2 J d d 



+ 8mk B T{d qi dq 2 ) 7(<?1 ^ U 9l dq 2 



Pt(qi,q2) (59) 



Since this equation is second order in q±,q 2 and their derivatives, we can propose a 
Gaussian ansatz, pt(qi, q-z) = exp{— c\q\ — c 2 q\q<i — c^q\ — ic\q\ — ic^q 2 — Cq}, without 
using the characteristic function associated with the Wigner function. This will lead 
to a set of coupled differential equations as in (120]) which is in principle solvable, but 
the basic difference from (1201) is that these equations will be non-linear due to the 
double derivative terms, e.g. 

d 2 

0^Pt(qi, ^2) = [(-2cigi - c 2 q 2 - ic^) 2 - 2ci] p t {q u q 2 ) (60) 

Thus they will be much more cumbersome to solve. 

In short, any equation of the form ( |59|) that contains double derivatives leads 
to a non-linear system of differential equations for q when we employ a Gaussian 
ansatz. This is the basic reason for using the characteristic function associated with 
the Wigner function to solve the evolution problem, where we solve a linear system 
of differential equations for q. See appendix [D] for more discussions. 

Let us now demonstrate how (1151) arises. The following identity follows from the 
Baker-Campbell-Hausdorff formula and will be useful in our calculations: 

e \(kq+xp) = e ±kq e jrxp gife = e j-xp e jrkq g-|f _ 
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Equation (|T3|) is derived as 

pt(k,x) =ti (p t ei^+^) 

du' (u'\p t e^ kq+xp) \u') 



du' J dv (u'\p t \v)(v\e% kq e^ xp e^ kx \u') 
d u ' e M«'+f) f dv (u'\p t \v)(v\u' + k) 



duei xu (u--\p t \u+-) , (62) 
{!«')} and {\v)} representing the momentum basis, and the inversion formula is 
{Pi\pt\P2) = / dud{u — ){u — \pt\u + — - — ) 

Znn J II 

1 f , -± X (E±+P2) f j ± xu/ P2-P1, 1 , P2-P1, 
1 dx e * y 2 J / due* xu (u — + — — ) 



<ix e" s3:1 ~ ' p t (p2 —pi,x) . (63) 



2nh 

The definition and the inversion formula for the position basis are derived similarly. 

To obtain f|T5|) from (TJ|), we will need the following identities which hold for opera- 
tors O and V and which follow from the basic commutation relations and the cyclicity 
of the trace: 

tr ([V, 0]e^ kq+xp) ^j = tr (o[e^ kq+xp \ V]) (64) 
[/(V), O] = cf'(V) if [V, 0\ = c a constant . (65) 
To calculate the first term on the right hand side of (|T5l) . note that 

9 ±ll.„_i_~ m \ d ( iu„ i„ m ikx\ i <Cl,_i_ m \ , k. 



— ( e f^P e * kq e-* 1 ?) = -Op - -)e^ {kq+xp) 
dx V J h 2 



which in turn implies 

d 



tr ({p, e>}e^ + ^)) = -2ih^-tx (Oe^ kq+xp ^j (67) 
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for any x- independent operator O. This gives 

tr ([^p 2 , p t ]e^A = _L tr (p,{p, [e^+^,p]} 



m ox V J m ox 



For the last term on the right hand side of (1151) . 



tr 



[q, {p, p t }]e^ + ^) = tr (p t {p, [ e *<**+-*>, g]}^ 
= x tr (p*{p, e^ fe9+a:p )}) 

= -2ih x— tr f p t e^ fe<z+ ^ ) = -2z^ x—p t (k, x) . (69) 
ax V / ax 

The second and the third terms are calculated in a similar way: 

tr ([g, [g, p t ]]e^+^) = tr (p,[[e^+^, g], g]) = x 2 p t (k,x) (70) 
tr ([p, [p,p 4 ]]e^ + ^) = tr (p t [[e^ + ^, p], p]) = fc 2 p 4 (A;,x) . (71) 

Putting all the terms together, <HM is obtained. 



B Green's Function for the Harmonic Oscillator 
System 

First, we calculate the contribution of the harmonic potential term to the differential 
equation for p t (k,x) in a way similar to ( 1681) : 

tr 



[\\muj 2 q \p t ]e^ +x ^ = l -mu 2 tr (p t {g, [ e *<***>, g]}^ 

= l -muj 2 x tr [p t {q,e^ kq+xp) ]^ 

= —ihmuj 2 x— tr ( p t e^^ kq+xp n = —ihmuj 2 x—p t (k, 



(72) 
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Together with this term, the equations for the coefficients in our Gaussian ansatz 
become: 



ci(t) 
c 2 (t) 
c 3 (t) 
c 4 (t) 



+ 



7 



m ' 8mkBT 

^^-2 1 c 2 {t)-2muj 2 c l {t) 
m 

2^mk B T 



h 2 
c 5 (t) 



47c 3 (t) - muj c 2 (t) 



(73) 



m 



= — 27c 5 (t) — muj 2 C4,(i) 
c 6 (t) = 

This system of equations is still linear but much more cumbersome to solve by hand. 
Hence, we used Mathematica® to solve it for the initial value 



1 



_ e -\(k-k ? e -i(2-z ) 2 



(74) 



We will not give the details of the calculation, which are quite lengthy. Upon rear- 
ranging the terms into complete squares and defining [i = a/7 2 — a; 2 , A t = 1 — e~ 2flt , 
we reach 



~Pt 0X0 ' ev (k,x) 



1 "7 

e 



ko—'- 



-(7-M)i 



At I 



x e 



!()- 



2^ ( (a* coth /it+7)fc— mu 2 i) 

"'' A i((^cothM-7)^+^) 



2 M 



- — f n Mi(t)fe7 



■M 2 (t)kx+m 2 uj 2 M3(t)x : 



(75) 



where 

Mi(i) = -- ^ \{e- 2 ^ cosh2 / ut - 1) 7 2 + r t u 2 + e" 27 * sinh2/it 7 /il 

[4(e- 27 * cosh 2/it - 1) 7 4 - 3(e^ 27 ' cosh 2/it - 1) 7 V 



16k 2 B T 2 fi 2 



+ T t ou 4 + Ae' 2lt sinh2/it 7 3 /i - e~ 27 * sinh2/it 7 w 2 /i ] 



M 2 (t) 



7 2 ( 7 -M)t A 2 

2/i2 6 



_^ 2 7 2 _ 
16k 2 B T 2 ^ 



+ 



[2(e~ 27 * cosh2 / ut - 1) 7 2 - (e~ 27 * + e" 27 * cosh2/it - 2) to 2 + 2e" 27i sinh2/it 7/1] 

(76) 
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MM) = At H e " 27 * cosh 2^ - 1) 7 2 - r t uj 2 + e~ 27t sinh2/it 7 J 

p.2 2 

- 16k 2 U T2fl 2 [( e " 27 * cosh2 ^ " !) 7 2 + r t uj 2 + e- 2 ^ sinh2^ 7/ x] . 

(77) 

Taking the limits e — > and rj — > 0, we obtain (j45]). 

C The Free particle as the Limit of the Harmonic 
Oscillator 

Note that the harmonic oscillator hamiltonian gives the free particle hamiltonian 
in the to — > limit, so we expect the Green's function and the density matrix of 
the harmonic oscillator to converge to those of the free particle in this limit. Thus, 
one can obtain the results for the free particle by first solving the problem for the 
harmonic oscillator and then taking the said limit. In our strategy, we will rather use 
this correspondence as an independent check of our results. Using Taylor series, we 
can expand Mi(t) around uj = 0: 

= MfrW + 4)-(3 -4e-^ + e-^) ^ 
4 7 2 

M 2 {t)= l -Y 2 t +0{u 2 ) (78) 
M 3 (t) = (l-e-^) + 0(uj 2 ). 

Note also that /i — ► 7 in the vanishing u limit. Then, it is trivial to recover (j25l) by 
inserting the above expressions into (|45|) and taking the uj — > limit. 



D The most General Case for the Gaussian Ansatz 



In this appendix, we discuss the most general equation for which the Gaussian ansatz 
can be employed. 
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Gaussian ansatz is applicable to any equation of the form 
df t (k,x) _ I" 

+ B k + C x 

+ dk + dx 

+ F k 2 + G kx + H x 2 

d d d d 

+ Lk— + M k— + N x— + P x— 
ok ox ok ox 

Q2 Q2 Q2 

+ Q — + R — — + S 



dk 2 dkdx dx 2 



ft(k,x) (79) 



We argued in appendix |A] that we have to solve a nonlinear system of differential 
equation unless Q,R, S = 0. One special case we can avoid nonlinearity is when 
F,G,H = 0. In that case, we can Fourier transform / in both k and x, and since the 
Fourier transform converts differentiation into multiplication, we do not have second 
order derivatives in the transformed equation. A single Fourier transformation can 
also be useful when F, R, S = and Q ^ 0, or H, Q, R = and S ^ 0. Roy and 
Venugopalan successfully use this approach in [9] to solve the time evolution of the 
harmonic oscillator density matrix for the non-Lindblad master equation, after certain 
change of coordinates in the position representation. However, if the "minimally 
invasive term" is introduced (which they do not attempt to do), their method cannot 
avoid having a second order derivative. For the rest of our discussion, we will set 
Q, R, S = and use the shorthand notation 

dft( £ x) = V(A, B, C, D, E, F, G, H, L, M, N, P)f t (k, x) (80) 

When we propose a Gaussian ansatz of the form (fTBl . we reach the following 
system of coupled linear equations 



c 2 (t) 
ca(t) 
c 4 (t) 
c 5 (t) 
V c e (t) J 



2L 
2N 


-2iD 





M 
L + P 

N 
-iE 
-iD 





2M 
2P 



-2%E 







L 

N 



\ 





M 

P 



iD iE / 



( c x \ 

c 2 
c 3 
c 4 
c 5 
V c 6 / 



+ 



/ ~F \ 

-G 
-H 
iB 
iC 
\ 'A ) 
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This is a system of inhomogeneous ordinary linear differential equations which 
can be solved by basic methods, but the dimension of the matrix makes the solu- 
tion intractable from a calculational point of view, even for mathematical software 
packages. 

The first observation that makes the calculation considerably easier is that Ci, C2, C3 
form an independent system of equations. This means, we can first solve for these 
three, then insert the solutions into the equations for C4 and C5 and solve the in- 
homogeneous equations for these two variables. We can finally insert C4, C5 into the 
equation for cq and find the solution by simple integration. This approach is tractable 
for Mathematica®, but the solutions are rather lengthy and give us little insight. 

The crucial step that simplifies (I8T1) is that by an affine transformation of the 
variables k, x in ([79]) . we can set the coupling terms D, E, M, N to for most cases, 
and have a diagonal matrix in ( IHTl) . Let us define the variables I, y such that 



which together with the scaling and swapping {k <-> x) can account for all linear 
transformations. This leads to the equation 



k = I + ay 
x = bl + y , 



(82) 



dfr(i,y) 



V(A', B>, C, D>, E>, F', G', H\ L', M', N', P')f t m (l, v) (83) 



dt 



with 



M' 



1 



(-6L + M - b 2 N + bP) 



N' 



1 — ab 
1 



(aL - a 2 M + N - aP) . 



(84) 



1 — ab 



By choosing 



a = 



{L-P) + y/(L - P) 2 + AMN 
2M 



b 



— {L - P) - y/(L - P) 2 + AMN 
2N 



(85) 



if M and N are both nonzero, and 



Approach to Thermal Equilibrium 



25 



if iV = (the case of M = is similar), we can set M' and N' to 0. The signs of 
the roots of the quadratics are chosen such that ab ^ 1, which ensures the linear 
independence of I and y. Note that this procedure cannot be used if MN = and 
L = P. 

Once we set M', N' = 0, given that V and P' are nonzero, we can shift our 
variables as 

, D ' 

m = t H 

L> 

z = y + §j, (87) 

which puts our equation into the form 

df ^ Z) ^ Z) = V{A", B\ C", 0, 0, F', G', H\ L', 0, 0, P')ti mz \m, z) . (88) 

This equation leads to six inhomogeneous ordinary differential equations which 
are not coupled, and thus can be solved quite easily. One can further simplify the 
equations if F' ^ 0, by scaling k \J~F'k to set the coefficient of the k 2 term to 1. 
By defining the function f[ mz ^ = f^^e^' 1 , the constant term A" can also be set to 
0. 

The special cases we did not discuss, e.g. MN = and P = L, can also be 
handled using similar techniques. Above transformations do not work when certain 
coefficients vanish or are equal to each other in (1791) . e.g N = and P = L. In 
these cases, solving flHTl) is already much easier before any affine transformation of 
the arguments of /. In short, using a Gaussian ansatz allows us to solve any equation 
in the form of ( 1791) without much trouble, as long as no nonlinear terms arise. 

We will not discuss the transformations of the coefficients in (!79|) under the affine 
transformations of the arguments of /, e.g C — ► C — > C" . They follow from basic 
algebra. 
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